Mesoscopic Shelving Readout of Superconducting Qubits in Circuit QED 
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We present a method for measuring the internal state of a superconducting qubit inside an on-chip 
microwave resonator. We show that one qubit state can be associated with the generation of an 
increasingly large cavity coherent field, while the other remains associated with the vacuum. By 
measuring the outgoing resonator field with conventional devices, an efficient single-shot QND-like 
qubit readout can be achieved, enabling a high-fidelity measurement in the spirit of the electron- 
shelving technique for trapped ions. We expect that the proposed ideas can be adapted to different 
superconducting qubit designs and contribute to the further improvement of qubit readout fidelity. 
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Superconducting nanocircuits [H, Q are considered 
promising candidates for diverse implementations of 
quantum information tasks [3]. In this context, circuit 
quantum electrodynamics (QED) [^,^], which studies su- 
perconducting qubits 0] coupled to on-chip microwave 
resonators, occupies a central role. To achieve the de- 
sired goals, it is important to implement high-fidelity 
two-qubit gates and efficient schemes to read out the 
qubit state In both cases, trapped- ion systems rep- 
resent the state-of-the-art for qubit realizations In 
particular, electron-shelving qubit readout has produced 
fidelity benchmarks of approximately 99.99% [105 . These 
astonishing achievements suggest the potential impact 
of transferring key ideas from quantum optics to circuit 
QED. Unfortunately, electron shelving relies strongly on 
the use of single-photon detectors [9], which are unavail- 
able in microwave technology in the range 1 — 10 GHz [ll| . 
Nevertheless, in this manuscript we show that a single- 
shot QND-like fast qubit readout can be designed by ex- 
ploiting the electron-shelving concept in circuit QED. 

We first present the physics of electron shelving in 
trapped ions. In Fig.[Tl we show a three- level atom where 
an unknown qubit state = a\g) + /^|e) is encoded in 
states |g) and |e). Via a laser beam, the ground state |g) 
is coupled to a third level |u), which can decay producing 
a continuous cyclic transition. In this case, the qubit is 
projected onto state |g) and many photons are emitted in 
free space, one at each cycle. In contrast, when the qubit 
is projected onto state |e), no photons are emitted. A 
lens is used to collect the photons more efficiently by im- 
proving the solid angle. Although the photodetector has 
a low efficiency 77d, the qubit readout fidelity can be very 
high. Typically, it is estimated through F = l — e~^^^, 
which rapidly approaches unity for r^d^ ^ 1, ^ being the 
number of detected photons. Key elements for electron 
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FIG. 1: (Color online) Sketch of electron shelving in trapped 
ions. The |g) ^ |u) transition is driven with a laser beam, per- 
forming a cyclic transition and emitting many photons when 
|g) is projected. No photons are detected when |e) is mea- 
sured. Undesired transitions are inhibited via selection rules. 



shelving are the use of three- level qubits, cyclic transi- 
tions, selection rules, and photodetectors. 

In the following, we present a method for implement- 
ing a single-shot QND-like fast high-fidelity readout of 
superconducting qubits. It preserves the spirit of elec- 
tron shelving, but it is suitably adapted to existent mi- 
crowave technology in circuit QED where, for example, 
single-photon detectors are unavailable. We assume that 
the qubit is prepared in an unknown pure state and that 
our task is to measure the spin operator . We consider 
a three-level superconducting qubit 12|, l3] inside an on- 
chip microwave resonator (acting as a cavity), as shown 
in Fig. O The initial qubit state is encoded in the two 
lower energy levels, |'0(O)) = a|g) +/^|e). In addition, we 
consider an anharmonic three-level qubit where the tran- 
sition frequencies are different, uOg^ ^ uo^^. Levels |e) and 
|u) are coupled resonantly to a resonator mode, but there 
is no dynamics because the resonator is initially empty 
and level |u) unpopulated. To start with the readout 
process, we drive the transition between levels |e) and 
|u) with a coherent resonant field with angular frequency 
c^d and amplitude ji transversal to the resonator axis. 

The system Hamiltonian, in the energy eigenbasis and 
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FIG. 2: (Color online) Sketch of the mesoscopic shelving qubit 
readout, (a) A three- level superconducting qubit is capaci- 
tively coupled (Cg) to a coplanar wave-guide microwave res- 
onator with angular frequency cjr and input and output ca- 
pacitors Cin and Cout, respectively. The qubit is also coupled 
to an orthogonal transmission line via . (b) The transition 
cjeu is resonant to the cavity and is driven with a transversal 
coherent field (magenta Une). Transition rates |e) |g) and 
|u} \g) are reduced by Purcell effect (green dashed lines) 
or selection rules (magenta dashed lines). 



after a rotating-wave approximation, can be written as 



H 



+nA (a^e-^^^^ + ae 



(1) 



Here, cr|^ = |u)(u| - |e)(e|, cr+ = |u)(e|, a^^ = |e)(u|, 
a (a^) are the bosonic annihilation (creation) operators 
of the resonator field, and ^eu, ^eu, and A are coupHng 
strengths. A describes the crosstalk between the driving 
field and the resonator, and its origin typically depends 
on the specific setup [14]. The qubit readout happens 
under the resonant condition c^r = = ^eu- We assume 
that the transition |u) ^ |e) is sufficiently long-lived such 
that it does not decay during the short operation time. 
Finally, our model considers enough energy anharmonic- 
ity so that the radiative decay rates associated with the 
transitions |e) |g) and |u) |g) are reduced by Purcell 
effect. Also, these transitions can be reduced exploiting 
the characteristic selection rules and symmetry breaking 
properties of superconducting qubits [13, J^]. 

We rewrite the Hamiltonian in a reference frame rotat- 
ing with the driving field frequency via the transforma- 
tion /7^°* = exp [—iuJd{ci^CL + o'^u^eu)'^] 5 obtaining 

H'""^ = ni]eu<u + %eu (^eti« + cr" a^) ^fiA (a^ + a) , (2) 

with (jg^ = cr+j + . We now apply the transformation 
= exp [—iftQuO'^uA under the strong-driving condition 
^eu ^ ^eu derive the effective Hamiltonian 
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-Kti + ^e"u)(« + + (at + a) . (3) 



The first part of the Hamiltonian simultaneously real- 
izes Jaynes-Cummings and anti-Jaynes-Cummings res- 
onant interactions. It does not generate Rabi oscil- 
lations, but conditional field displacements [17], while 
the second term implements a resonant displacement. 
The initial qubit-field state is |'0(O)) = a|g)|0) + 
+ I-)) with a^J±) = ±\±). After an in- 

teraction time t, the state is 

\m) = a\gMt)) + ^[\+Mt) + m) 

+\-)\-m + m)]- (4) 

Here, the coherent states | ± f]{t) + ^(^)), with 
f]{t) = —igQ^t/2 and v{t) = — zAt, are generated 
by the displacement operators V[±.fi{t) + ^if)] = 
exp {[±f]{t) + v(t)\a^ - [±f]*(t) + v*(t)\a). In general, 
we expect the crosstalk to be small, so that A <C ^eu/2 
and V{t) <C f]{t). When the measurement starts, the 
applied driving field yields many intracavity photons 
^fn(^) ^ with probability \a\^ if the state |e) is 

projected. If the state |g) is selected, it yields a few pho- 
tons Nf^{t) = \iy{t)\^ < \f]{t)\^ with probability \f3\^. 

We now add to our model a zero-temperature dissi- 
pative reservoir for the cavity field, characterized by a 
decay rate k,. The corresponding master equation reads 



Pq-f = -■^[^efr,Pq-f] + A Pq-f 



(5) 



with >CfPq_f = A>:L[a]pq_f such that 

>^fPq-f = (2apc^-fa^ — a^apc^-f — pc^-fa^a) /2 (6) 

and expansion pq_f(t) = \j){k\^pf (t). Here, 

it is possible to find analytical solutions for pj^j t) = 
{j\pq-i{t)\k) using standard phase-space tools [20| and 
the method of characteristics to solve the partial differ- 
ential equations [2l|. The solutions read 
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FIG. 3: (Color online) Intracavity mean photon number for 
the mesoscopic shelving readout of a CPB in the charge-phase 
regime with conservative parameter set: Ec/h — Ej/h = 

27T X 10 GHz, ^eu = ^eu/S = 27T X 150 MHz, K = 27TX 1.6 MHz, 

= 107ge = 2 MHz, 7gu = 0, and = 2 MHz. The dotted 
and dashed curves correspond to the analytical results and the 
solid lines to the numerical results. Note that, in the absence 
of losses, there is still a difference between the analytical and 
the numerical results. This is due to the off-resonant couplings 
and multilevel character of the realistic model. 

with r]{t) = -igen/i^[l-e-^'^^], v(t) 
— 2iA//^ [l — e~'^^/^] . For a small crosstalk A, the 
leakage rate of outgoing photons N^^^(t) when state |e) 
is measured can be estimated as 

where N^^{t) = \r]{t)\'^ is the intracavity mean photon 
number. A^out(^) gi'ows very fast well below decoherence 
times. It can be measured, e.g., by means of a data 
acquisition card, which follows a phase-preserving or even 
a more quiet phase-sensitive jlSj linear amplifier. We 
also notice that one can profit from the generated large 
intracavity field to adapt to other readout techniques [l9| . 

The physical concepts behind the mesoscopic shelving 
are general and can be adapted to different qubits and se- 
tups. We exemplify here with a possible adaptation to a 
Cooper-pair box (CPB) coupled to a microwave resonator 
of angular frequency ujr- Here, the CPB has a Josephson 
energy Ej and charging energy = (2e)^/2Ctot, where 
Ctot is the total island capacitance. We refer to a system 
that is essentially the one in Ref. [3], with the addition 
of a transmission line, orthogonal to the resonator, for 
driving the qubit (see Fig. [2]). Using n, the number op- 
erator for excess Cooper pairs on the CPB island, and (p, 
the phase difference across the Josephson junction, the 
Hamiltonian can be written as 

H = Ec{n — Ux)'^ — Ej cos (p -\- hcOrCt^ CL , (10) 



where 

(2e)n, = + Cl'V^^^t) + qFo(a^ + a) . (11) 

Here, Cg (/ = {dc,ac,r}) are effective gate capacitances, 
V^'^ is the gate voltage that defines the working point (we 
choose the so-called "sweet spot" C^^l/g^72e = 1/2), V^'^ 
is the voltage of the orthogonal driving field, a {a^) refers 
to the cavity field, and Vq is the resonator zero-point volt- 
age. Note that the CPB is coupled to the resonator and 
the ac drive via the charge number operator. The clas- 
sical gate charge C^^V^^ and the quantum gate charge 
CgVb represent small deviations from the sweet spot. 

We can rewrite in a basis of CPB eigenstates re- 
stricted to the first three energy levels, the ground state 
|g) and the first and second excited states, |e) and |u), 
respectively. This leads to an effective Hamiltonian for 
the driven qubit-resonator system, 

H = Ho -\- Hint + Hd , (12) 
Ho = Ej\j){j\^hw,a^a , 

J=g,e,u 

Hint = f^(^gu|u)(g| +^eu|u)(e|)(a''' +a) + H.c. , 
Hd = (coscjdt) ^ {h^jk\j){k\ + H.c.) . 

j</e=g,e,u 

The coupling strengths gjk = {Ec/he)njkC^Vo and 
rtjk = {Ec / he)njkC^^V^^ (t) are proportional to the ma- 
trix elements njk = {k\n\j). In order to obtain the time 
evolution of the complete system, including the relax- 
ation and dephasing of qubit transitions, we numerically 
solve the master equation for the qubit-resonator density 
matrix p{t) 

P = -^[H,p]^ Cfp^ CdP ^ (13) 
where, using the functional L defined in Eq. ([5]), we have 

c^p = j2i,,L[\k){j\]p+j2:^L[\j){j\]p. (u) 

j<k j 

Here, hi is the decay rate of the resonator, jjk {j^k = 
{g, e, u}) are the relaxation rates for the transitions |e) 
|u) and |u) |g), and is the dephasing rate, which 
we take to be equal for all coherences. In the numerical 
solution of Eq. ([13]), we truncate the resonator Hilbert 
space to 25 photon number states due to technical lim- 
itations. In addition, we make sure that the population 
of the fourth qubit eigenstate is negligible. Clearly, the 
condition 7eu + 7gu <^ ^eu is crucial for our method to be 
efficient. Therefore, a qubit layout with suppressed 7gu 
is preferred for the shelving readout. 

The results for the intracavity mean photon number 
with a conservative set of parameters, in the analytical 
and numerical cases, are shown in Fig. [3l We see that. 
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although the fuh dynamics in Eq. ([T3|) is considerably 
more complex than the one in Eq. ([5]) , the simple analyt- 
ical model captures the essence of the system dynamics. 
The main influence of a realistic description is a small 
reduction in the intracavity mean photon number. We 
observe that, given the short interaction times displayed 
in Fig. [3l the resonator decay rate Hi alone has a small 
effect in the cavity population, while the finite lifetime 
of states |e) and |u) is slightly more important. In this 
manner, we feel comfortable to extrapolate the analytical 
results for the cavity population including cavity losses 
for short measurement times to make further estimations. 

The signal-to-noise ratio (SNR) after a measurement 
time Tm is the ratio between the accumulated num- 
ber of outgoing photons and the accumulated noise 
The latter is dominated by the amplifier noise, namp = 
k^T^/hw^ — 25, where is its associated noise temper- 
ature. In this manner, 

SNR(..) = (15) 

where B = max{/^, 7j/c} is the measurement bandwidth. 
We now estimate the SNR for three relevant consecu- 
tive times. First, we use the maximum simulated time 
8 ns, corresponding to 10 intracavity photons (cf. 
Fig. E]). We obtain a SNR 0.2. Considering that 
our simulations include all relevant system details with- 
out any approximation ^2[, this is a remarkable result 
for such an extremely short measurement time. Using 
our analytical results including resonator dissipation, see 
Eq. (|9]), we estimate that a critical measurement time 
r^^^ ^ 19 ns is necessary to reach the condition SNR ~ 1. 
This is the minimum time required for a single-shot mea- 
surement of the qubit state |e). Finally, to achieve high- 
fidelity qubit readout, we choose the measurement time 
:^ 50 ns which corresponds to SNR 6.2 and fidelity 
F = 1 - e-^^^ = 0.998. Notably, <C l/7ge. Con- 
sequently, we expect a single-shot measurement of the 
qubit state |e) with fidelities close to 1. The proposed 
mesoscopic shelving qubit readout is of a QND-like char- 
acter, due to the continuous cavity field amplification 
in each measurement event. In addition, is at least 
one order of magnitude shorter than typical measurement 
times employed in the state-of-the-art experiments based 
on dispersive readouts. Note that, even for a driving- 
resonator crosstalk of A = 27r x 10 MHz [3], the cavity 
population associated with the measurement of state |g) 
is well below the amplifier background noise level. 

In summary, we have presented a novel qubit readout 
scheme based on a mesoscopic shelving technique, allow- 
ing a fast high-fidelity single-shot QND-like measurement 
of superconducting qubits in circuit QED. 
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